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; Abstract 

We have extracted the constituent contributions to the spin of the proton 
from recent data at CERN, SLAC and DESY. The valence, sea quark and 
antiquark spin-weighted distributions are determined separately. The data 
appear to imply a small to moderate polarized gluon distribution, so that 
\ the anomaly term is not significant in determining these contributions. We 

have analyzed the consistency of the results obtained from various sets of 
data and the Bjorken Sum Rule. All data are consistent with the sum 
rule, but they differ in the contribution of the strange sea to proton spin. 
This and the remaining uncertainty in the polarized gluon distribution pose 
' unanswered questions about hadronic spin. Only further experiments which 

extract information about the polarized gluon and sea will reconcile these 
differences. We suggest specific experiments which can be performed to 
determine the size of the polarized sea and gluons. 
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CL,; Introduction 

Oh! 

One of the goals of high energy spin physics is to determine the contribu- 
tions of quarks and gluons, as well as the effect of the orbital motion, to nucleon 
spin. Significant interest in high energy polarization was generated when the 
European Muon Collaboration (EMC) 1 analyzed polarized deep-inelastic lepton- 
hadron scattering (DIS) data, which implied that the Bjorken sum rule (BSR) 
of QCD 2 was satisfied and the Ellis-Jaffe sum rule 3 based on a simple quark 
model was violated. Recently, the Spin Muon Collaboration (SMC) group from 
CERN 4 , the experimental groups from SLAC 5 and the HERMES group at DESY 6 
measured the polarized structure functions g\ to low x and have added the cor- 
responding neutron and deuteron structure functions and gf. They also im- 
proved statistics and lowered the systematic errors from the EMC data. 

An advantage to using spin in probing hadronic structure is that the theo- 
retically calculated spin-weighted parton distributions are related to the experi- 
mentally measured cross sections by the polarized structure functions, g\. The 
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measured asymmetry in deep-inelastic lepton-hadron scattering (DIS) is given 
by: 
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= + ijA 2 ). 
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The spins of the target and beam are aligned or anti-aligned with the momentum 
of the infraction in measuring these cross sections. The terms D and rj are known 
kinematic factors for each experiment. The structure function g± is normally 
extracted from the asymmetry A\ by using the approximation 

2 _ ^(x)F 2 (x,Q 2 ) 
5l(X ' Q) 2x(l + J R) ' (2) 

where F2 is the corresponding unpolarized singlet structure function and R = 
\J{<jl/vt), the ratio of the cross sections for absorbing longitudinal and transverse 
virtual photons. It is assumed that the transverse part of the asymmetry A\ is 
small and that A\ is relatively independent of Q 2 , which has been implied by 
these experiments. 

The integrated polarized structure function, JP( n ) = Jq g^ n \x) dx, is related 
to the polarized quark distributions by 

pin) = J_ (1 _ a «Tr)<[4(i) Aut) + 1(4)A4 + 4(l)(Au a + An) + l(4)(Ad a + Ad) 

lo 

+ (As + As)]). (3) 

The QCD corrections, characterized by a c ° rr , have been calculated to O(o^) 7 ' 8 
and are 

a corr ^ ^ + 3. 58 33(^£)2 + 20 .2153(— ) 3 + 130(— ) 4 , (4) 
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where the last term is estimated. Equations (1) through (4) provide a direct 
means to extract information about the polarized quark distributions from the 
DIS experiments. There have been some recent theoretical approaches to this 
problem. 7,9 We have done a detailed flavor dependent analysis including the QCD 
corrections and the effect of the gluon anomaly. It is assumed that the polarized 
gluon distribution is of small to moderate size and we determine the resulting 
polarized quark distributions for each set of data using the appropriate sum 
rules. The key elements of our approach are: 



• determine the valence contribution to the spin using the BSR 



• find sea integrated parton distributions for each flavor by breaking the 
SU(6) symmetry with the strange quarks and using the sum rules with 
data as input 

• include higher order QCD corrections and the gluon anomaly for each flavor 

• discuss similarities and differences between the phenomenological implica- 
tions of the different experimental results, and 

• suggesting a set of experiments which would distinguish the quark and gluon 
contributions to the proton spin. 

This approach differs from that of others in that we use sum rules in con- 
junction with a single experimental result to extract the spin information and we 
break the flavor symmetric sea while including anomaly contributions. 



Valence Quarks 

Fundamentally, we assume that the proton is comprised of valence quarks, 
whose integrated polarized distribution is given by: (Aq v (Q 2 )}. We construct 
the polarized valence quark distributions from the unpolarized ones by starting 
with a 3-quark model based on an SU(6) proton wave function. The valence 
quark distributions can be written as: 



where cos#d is a "spin dilution" factor which vanishes as x — > and becomes 
unity as x — > 1, characterizing the valence quark helicity contribution to the 
proton. 10,11 The spin dilution factor is adjusted to satisfy the Bjorken Sum Rule 
(BSR). The BSR relates the polarized structure function gi(x), measured in po- 
larized deep-inelastic scattering, to the axial vector current A3, which is measured 
in neutron beta decay. This sum rule is considered to be a fundamental test of 
QCD. In terms of the polarized distributions and our assumptions about the 
flavor symmetry of the u and d polarized sea, the BSR can be reduced to: 
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Thus, the valence contributions are determined uniquely by this model. The va- 
lence distributions are not sensitive to the unpolarized distributions used to gen- 
erate them. 12 ' 13 With our values (Au v ) = 1.00 ± 0.01 and (Ad v ) = -.26 ± 0.01, 
both the BSR and magnetic moment ratio, fip//J, n ~ —3/2 are satisfied. This 
results in a spin contribution from the valence quarks equal to 0.74 ± 0.02. The 
quoted errors arise from data on A% = gA/gv, and the differences in choice of the 
unpolarized distributions. 

Sea Quarks 

The proton is also filled with a quark sea, whose lightest flavors should dom- 
inate the spin, since the heavier quarks would be significantly harder to polarize. 
We assume that the quark and antiquark flavors are symmetric, but break the 
SU(6) symmetry of the sea by assuming that the polarization of the heavier 
strange quarks is suppressed. 10 The sea distributions are then related by: 

Au(x, Q 2 ) = Au(x, Q 2 ) = Ad(x, Q 2 ) = Ad(x, Q 2 ) 

= [1 + e] As(x, Q 2 ) = [1 + e]As(x, Q 2 ). (7) 

The e factor is a measure of the increased difficulty in polarizing the strange 
quarks. The DIS data are used to determine e and the overall size of the polarized 
sea. Additional constraints are provided by the axial-vector current operators, 
^8 an d Aq. 

The coefficient As is determined by hyperon decay and is related to the po- 
larized quark distributions by: 

A 8 = {[Au v + Ad v + Au s + Au + Ad s + Ad - 2As - 2 As]) « 0.58 ± 0.02. (8) 

Aq is related to the total spin carried by the quarks in the proton. We can relate 
these axial currents and the structure function g\ in the anomaly-independent 
form: 

A = 9(1 - a?"')" 1 C &(x) dx - -A % - -A z « (Aq tot ). (9) 
J0 4 4 

Gluons 

The gluons are polarized through Bremsstrahlung from the quarks. The inte- 
grated polarized gluon distribution is written as: (AG) = Jq 1 AG(x,Q 2 ) dx. We 
cannot determine a priori the size of the polarized gluon distribution at a given 
Q 2 . The evolution equations for the polarized distributions indicate that the 
polarized gluon distribution increases with Q 2 and that its evolution is directly 
related to the behavior of the orbital angular momentum. 14 Thus, one assumes a 



particular form for the polarized gluon distribution at a given Q 2 and checks its 
consistency with data which are sensitive to AG(x,Q 2 ) at a particular Qq. 

The model of AG that is used has a direct effect on the measured value of 
the quark distributions through the gluon axial anomaly, 15 which has the form: 

r(Q2) = N f Os(Q 2 ) £ AG(Xj Q2) (1Q) 

where Nf is the number of quark flavors. For each quark flavor, the measured 
polarization distribution is modified by a factor: (Aqi) — T(Q 2 )/Nf. Thus, the 
quark spin contributions depend indirectly on the polarized gluon distribution. 
In a naive quark model, (Aq) = 1 and AG may be quite large to be consistent 
with data. 16 However, a reasonably sized AG is possible if the sea has a suitably 
negative polarization. We have considered two possible models for calculating 
the anomaly: (1) AG = xG (indicating that the spin carried by gluon is equal 
to its momentum) and (2) AG = 0, which is an extreme case for bounding 
the distributions. We believe that the present data imply that anomaly effects, 
and thus the overall integrated polarized gluon distribution, is limited at these 
energies. 

The polarized distributions are related to the orbital angular momentum of 
the constituents by the J z = \ sum rule: 

J z = \ = l -(Aq v ) + \(AS) + (AG) + L z . (11) 

The right hand side represents the decomposition of the constituent spins along 
with their relative angular momentum, L z . Although this does not provide a 
strict constraint on either Aqtot or AG, it does give an indication of the angular 
momentum component to proton spin. 



Phenomenology 



We use SMC 4 , SLAC 5 and DESY 6 data to extract information about the 
flavor dependence of the sea contributions to nucleon spin. We can write the 
integrals of the polarized structure functions, Jq 1 g\ dx in the terms of the axial- 
vector currents as: 



p= f'g^dx 
Jo 

= f gi(x)dx 
Jo 
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where lod is the probability that the deuteron will be in a D-state. Using N-N 
potential calculations, the value of ur> is about 0.058. 17 The BSR can then be 
used to extract an effective I p value from all data using the form of equation (12) 
above. 

Since the anomalous dimensions for the polarized distributions have an ad- 
ditional factor of x compared to the unpolarized case, early treatments of the 
spin distributions assumed a form of: Aq(x) = xq(x) for all flavors. 18 We have 
compared this form of the distributions to those extracted from the recent data, 
using the defined ratio ri = ^ qsea ^ ex p f or each flavor. Any deviation from r? = 1 
would indicate that the simple model for generating the polarized distributions 
is inaccurate. 

In order togenerate the x-dependent distributions, we have used the unpo- 
larized distributions 12 ' 13 with our extracted value of rj and the assumption that: 
Aq(x) = 7]xq{x) for each of the sea flavors. For the valence distributions, we 
have used equation (5) with the dilution factor of reference 10. There is no rea- 
son a priori to suspect that a global fit to the integrated distributions should 
imply a satisfactory x-dependent fit to the data. However, our results indicate 
that this form gives very good x-dependent parametrizations for the polarized 
distributions. 

The analysis (for each polarized gluon model) proceeds as follows: 

• Extract a value of I p from either the data directly or via the BSR in the 
form of equation (12), 

• use equation (9) to extract Aq. Then the overall contribution to the quark 
spin is found from (Aq tot ) = Aq + T. 

• Use the value As from the hyperon data with equations (8) and (9) to 
extract As for the strange sea, 

• find the total contribution from the sea from (Aq tot ) = (Aq v ) + (AS), 

• determine e and the distributions (Au) sea = langleAd) sea from equation 
(3) and the strange sea results. 

• Finally, the J z = 1/2 sum rule gives L z . 



Results for the integrated distributions are given in Table I. 



Table I: Integrated Polarized Distributions: 
AG = xG (above line), AG = (below line) 



Quantity 


SMC(P) 


SMC(I d ) 


E15A(I n ) 


EU3(I d ) 


HERMES 

(r) 


< Au >sea 


-.077 


-.089 


-.063 


-.068 


-.050 


< As > 


-.037 


-.048 


-.020 


-.028 


-.010 


< Au > tot 


0.85 


0.82 


0.87 


0.87 


0.90 


< Ad > tot 


-.42 


-.43 


-.39 


-.40 


-.36 


< As > tot 


-.07 


-.10 


-.04 


-.06 


-.02 


Vu = Vd 


-2.4 


-2.8 


-1.9 


-2.1 


-1.5 


Vs 


-2.0 


-3.0 


-1.2 


-1.6 


-0.6 


e 


1.09 


0.84 


2.10 


1.41 


4.00 


r 


0.06 


0.06 


0.08 


0.08 


0.07 


p 


0.136 


0.129 


0.134 


0.131 


0.135 


< Aq > tot 


0.36 


0.29 


0.45 


0.41 


0.52 


< AG > 


0.46 


0.46 


0.45 


0.44 


0.44 


L z 


-.14 


-.211 


-.18 


-.15 


-.22 



< Au >tot 


.83 


.80 


.85 


.84 


.88 


< Ad > tot 


-.44 


-.45 


-.41 


-.43 


-.39 


< As > tot 


-.09 


-.12 


-.07 


-.08 


-.04 


IP 


.136 


.129 


.134 


.131 


.135 


< Aq > tot 


0.30 


0.23 


0.37 


0.33 


0.45 


r 


0.00 


0.00 


0.00 


0.00 


0.00 




0.35 


0.39 


0.32 


0.35 


0.28 



From the results in Table I, it is obvious that the naive quark model is not 
sufficient to explain the proton's spin characteristics. Nor is the simple model 
for extracting the polarized distributions accurate. Thus, data have indicated 
that we must modify our initial assumptions regarding constituent contribution 
to proton spin. Some conclusions which can be drawn from the data are: 



(1) The total quark contribution to proton spin is between 1/4 and 1/2. The 
errors in generating these results are due mostly to experimental errors and deter- 
mination of which model of the polarized gluons to use. Thus, the uncertainties 
related to the quark spin content and the size of AG are comparable. 

(2) Considerable discussion regarding these measurements focuses on the Ellis- 
Jaffe sum rule (EJSR), 3 which predicts the values of g\ and g\ using an unpolar- 
ized strange sea. This differs considerably with the analysis of these data. The 



up and down sea contributions seem to agree within a few percent. However, all 
of the proton and deuteron data imply a larger polarized sea with the strange sea 
polarized greater than the positivity bound. 19 It is interesting to note that the 
results obtained from the SMC proton data are consistent with a recent lattice 
QCD calculation of these parameters. 20 The results among the different experi- 
ments preclude us from extracting specific contributions from each flavor to the 
proton spin with any degree of certainty. However, the results from the various 
data can be categorized into distinct models, characterized by the size of the 
non-zero polarized sea. 

(3) The values of r\ deviate considerably from unity for most of the data, im- 
plying that the relation between unpolarized and polarized distributions is likely 
more complex than originally thought. 

(4) This analysis implies that the anomaly correction is not large. If the 
anomaly term were larger, due to a large AG, the strange sea would be positively 
polarized, while the other flavors are negatively polarized. There is no known 
mechanism that would allow this cross polarization of different flavors. Thus, we 
conclude that these data imply that the polarized gluon distribution is of small to 
moderate size. The key conclusions from the zero and small anomaly models are 
not significantly different. Further, even if there are higher twist corrections to 
the anomaly at small Q 2 , the anomaly will not reconcile differences in the flavor 
dependence of the polarized sea. 21 

(5) The orbital angular momentum extracted from data is also much smaller 
than earlier values obtained from EMC data. 16 A point of interest is that the 
zero AG model implies a positive orbital angular momentum, while the other 
model gives a negative result. Thus, although L z is likely small, its sign is still 
in question. 

(6) The extracted value for I p is comparable for all data and well within 
the experimental uncertainties. This implies agreement about the validity of the 
Bjorken Sum Rule. We have arrived at this conclusion by using the BSR to ex- 
tract an effective I p , in contrast to the experimental groups, which used data to 
extract the BSR. There is general agreement that the BSR (and thus QCD) is in 
tact. 

Clearly, these experiments have contributed to the progress of understand- 
ing the relative contributions of the constituents to the proton spin. They have 
probed to smaller x values, while decreasing the statistical and systematic er- 



rors. This, coupled with theoretical progress in calculating higher order QCD 
and higher twist corrections have allowed us to narrow the range of these spin 
contributions. Although the flavor contributions to the proton spin cannot be 
extracted precisely, the range of possibilities has been substantially decreased. 
The main differences are the questions of the strange sea spin content and the 
size of the polarized gluon distribution. We stress that more experiments must be 
performed to determine the relative contributions from gluons and various flavors 
of the sea. 

The x-dependent distributions are in very good agreement with proton, neu- 
tron and deuteron data. The graphs which follow compare x binned data with 
the polarized x-dependent distributions generated from results on i] and equation 
(3). These were found using the GRV 12 and the MRS 13 unpolarized distributions. 
The differences between the two sets of distributions are at small-x, where the 
data is most uncertain. More DIS experiments should be performed at small-x 
to distinguish between models and to address the controversy regarding which 
contributions to g\ are dominate in this kinematic region. 

Possible Experiments 

There are a number of experiments which are technologically feasible that 
would supply some of the missing information about these distributions. Detailed 
summaries can be found in references 22 and 23. The large average luminosities 
of these experiments and the success of Siberian Snakes makes all of the following 
feasible. 

Deep Inelastic Scattering: The E155 experiment has been approved at SLAC. 
These experiments are designed to probe slightly smaller x while greatly im- 
proving statistics and systematical errors. With lower error bars at small x, the 
extrapolation should achieve a more accurate value for the integrated distribu- 
tions and narrow the ranges of constituent spin contributions even further. 

Lepton Pair Production (Drell-Yan): The Relativistic Heavy Ion Collider 
(RHIC) at Brookhaven is designed so that polarized pp and pp experiments can 
be performed at large energy and momentum transfer ranges. The energy range 
will be covered in discrete steps of about 60, 250 and 500 GeV, but the momentum 
transfer range covers 0.005 < Q 2 < 6.0 GeV 2 in a fairly continuous set of steps. 
The PHENIX detector is suitable for lepton detection and the wide range of 
energies and momentum transfers could yield a wealth of Drell-Yan data over a 
wide kinematic range. The x-dependence of the polarized sea distributions could 
then be extracted to a fair degree of accuracy. 

Jets, pions and direct photon production: The SPIN Collaboration proposes 
a set of experiments, which are in the kinematic region where the measurement 
of double spin asymmetries in jet production would give a sensitive test of the 
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Figure 1: The x- dependent structure function g\ is compared with data 
as a function of x. 




Figure 2: The z-dependent structure function g" is compared with data 
as a function of x. The solid line represents the GRV generated distri- 
butions and the dashed line the MRS generated distributions. 
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Figure 3: The x- dependent structure function gf is compared with data 
as a function of x. The solid line represents the GRV generated distri- 
butions and the dashed line the MRS generated distributions. 



polarized gluon distribution's size. Naturally, this measurement has an effect on 
both AG and the anomaly term appearing in the polarized quark distributions. 24 
The STAR detector at RHIC is suitable for inclusive reactions involving jet mea- 
surements, direct photon production 25 and pion production. All of these would 
provide excellent measurements of the Q 2 dependence of AG since all are sensi- 
tive to the polarized gluon density at differing Q 2 values. Charm production in 
polarized collisions are also sensitive to AG and should be performed at RHIC. 
Should DESY proceed with plans to polarize their proton beam, many of these 
experiments could be performed there, complementing the kinematic regions cov- 
ered by RHIC and CERN. 

There has been considerable discussion about performing the COMPASS po- 
larization experiments at the LHC at CERN. Depending on the approved exper- 
iments, there is the possibility of probing small x and doing polarized inclusive 
experiments to measure both sea and gluon contributions to proton spin. These 
could be made in complementary kinematic regions to those of the other acceler- 
ators. There are tentative plans to do polarized production, which provides 
a measure of the x-dependent sea distributions. Polarized H /± production is also 
planned at SLAC and would provide useful sea information in a slightly different 
kinematic region than that of CERN. 

Tests of the valence quark polarized distributions can be made, provided a 
suitable polarized antiproton beam of sufficient intensity could be developed. 26 
This would provide a good test of the Bjorken sum rule via measurement of (Aq v ) 
and the assumption of a flavor symmetric up and down sea. This should be an 
experimental priority for the spin community. 

Existing data indicate that the spin structure of nucleons is non-trivial and 
has led to the formulation of a crucial set of questions to be answered about 
this structure. The experiments discussed above can and should be performed in 
order to shed light on the spin structure of the nucleons. 
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